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Abstract. One of the ways that analysis on fractals is more complicated than 
analysis on manifolds is that the asymptotic behavior of the spectral counting 
function N{t) has a power law modulated by a nonconstant multiplicatively 
periodic function. Nevertheless, we show that for the Sierpinski gasket it is 
possible to write an exact formula, with no remainder term, valid for almost 
every t. This is a stronger result than is valid on manifolds. 



1. Introduction 
The Weyl asymptotic formula for the counting function 

(1.1) Nit) = #{A, < t} 

for the eigenvalues {Aj} of — A on a compact manifold fl (with or without boundary) 
states that 

(1.2) N{t) - c„vol(^2)r/2 as t ^ oo 

where c„ is a constant depending on the dimension n of fl. There are more precise 
statements such as 

(1.3) N{t) - cn voi(ri)r/2 + o(t("-i)/2) 

for manifolds without boundary or with smooth boundary, and 

(1.4) Nit) = c„ vol(f7)t"/2 + 0(i'^/2) 

for manifolds with non-smooth boundary of Minkowski dimension d [L] . 

In the case of Laplacians on fractals, the situation is more complicated. For 
example, for the Sierpinski gasket iSG) with either Dirichlet or Neumann boundary 
conditions, the analogous statement is 

(1.5) iV(i) - G(i)r as i -> cx) for a = 

logs 

where G is a multiplicatively periodic function 

(1.6) G(5i) = Git) 

that is bounded away from zero, and discontinuous (hence non-constant). This was 
first observed by Fukushima and Shima [FS] based on a precise description of the 
spectrum. It was later put into a wider context by Kigami and Lapidus [KL] based 
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on the Reneval Theorem. Using a refinement of the Reneval Theorem, Kigami [Kl] 
showed that the error is bounded, 

(1.7) N{t) = G{t)e + 0{l). 

Note that this is consistent with the fact that the boundary of SG consists of three 
points, hence has dimension zero. Similarly, the explanation for the value of a is 
that a = d/{d + 1) where d = log 3/ log 5 is the dimension of SG in the resistance 
metric and d + 1 is the order of the Laplacian [S2] . 

We already see that (1.7) is stronger than (1.3). The point of this note is that 
even more is true; we can have an exact formula with no remainder term at all, 
provided we restrict attention to almost every t. More precisely, let SG denote the 
double cover of SG (two copies glued together at corresponding boundary points). 
Then SG is an example of a fractafold without boundary [SI]. Neumann eigen- 
functions on SG extended evenly to SG and Dirichlet eigenfunctions on SG extend 
oddly to SG give a complete set of eigenfunctions on SG. In particular 

(1.8) N{t) = N^it) + Noit) 

where N, Nn and Nd denote the eigenvalue counting function on SG, and SG 
with Neumann or Dirichlet boundary conditions, respectively. 

Theorem: There exists an open set A c (0, oo) whose complement has measure 
zero, with 5 A = A, such that for each t £ A there exists mo{t) such that 

(1.9) iV(5"t) = 2G(5'"t)(5'"i)« = 2G(i)i«3'" 
for all m > mo{t), and G is continuous on A. Also 

(1.10) iVAr(5™i) = G(i)i"3" + Gi(i) 

7V^(5™i) = G(0t"3"-Gi(i) 

for all m > mo{t), where Gi{5t) = Gi{t) and G\ is hounded. 

Our proof is based on the "spectral decimation" description of the spectrum 
from [FS] . It seems likely that a similar result holds for many pcf fractals for which 
spectral decimation is valid ([Sh], [B-T]), but this would require a different proof. 
It is straightforward to extend the result to all fractafolds without boundary based 
on SG [SI]. 

It is instructive to compare our result with the case of the second derivative 
on an interval. To be specific, consider the interval [0, tt] with Dirichlet boundary 
conditions. The eigenvalues are the perfect squares and N{t) = [t^^"^]- Then we 
have N{t) = t^^"^ exactly when t is a perfect square, so N{4"^t) = 2™t^/^ for ah t of 
the form 4~"fc^ for m > n. Thus the analog of the set ^ is a countable dense set. 
This is not the first example of a result in fractal analysis that is stronger than the 
corresponding result in classical analysis (see [S3] for convergence of Fourier series, 
and also [BK], [BS], [S5]). We hope it will not be the last example. 

The reader is referred to the books [K2] and [S4] for any undefined notation. 

2. Structure of the Spectrum 

We review briefly the structure of the spectrum for SG and SG with Dirichlet 
or Neumann boundary conditions ([FS], [GRS], [SI], [S4]). The lowest eigenvalue 
is of course zero, corresponding to the constant eigenfunction, and this occurs with 
multiplying one for SG and SG with Neumann boundary conditions. After that. 



EXACT SPECTRAL ASYMPTOTICS ON THE SIERPINSKI GASKET 



3 



the eigenvalues arrange themselves in cycles Ci, C2, . . .. Each cycle contains ten of 
distinct eigenvalues. The eigenvalues and their multiplicities in a cycle are shown 
in Table 1. We factor out the powers of 2 to write uniquely 

(2.1) fc = 2^(2^-1) 

where j is a nonnegative integer and £ a positive integer. There are increasing 

families of primitive eigenvalues {\n^}, {^n^} and {aI^^} for n = 1,2,3, The 

upper index denotes the eigenvalue of the graph Laplacian of the restriction of the 
corresponding eigenfunctions to the nine element graph Fi of one approximating 
SG (Figure 2.1). 



Figure 2.1. The graph Fi (dotted lines indicate identified points) 



Table 1. Eigenvalues and their multiplicities in the cycle Ck 



Eigenvalues 


Mn 


Md 


M 


Nn-Nd 


'^2k-l 





1 


1 





A (3) 
'^2k~l 


2 





2 


2 


U5) 
'^2k-l 





2 


2 







3 





3 


3 


x(5) 
^2fe 





2 


2 


1 


'^2k 


2 





2 


3 


'^2k 





1 


1 


2 




.V + i-l 

2 


3J+H3 
2 


3^+1 + 1 







3-'+^+3 


3^+^-3 




3 


« + l \(5) 
A2^_i 


3^ + 1-1 
2 


3^ + ^+3 
2 


3^+1 + 1 


1 



The table lists the eigenvalues in increasing order (top to bottom), followed by 

their multiplicities in each of the three cases. Note that M = Mm + Mu- The last 
column gives the difference of the eigenvalue counting functions N]\[{t) — Njy{t) for 
t equal to the eigenvalue. Since this is a cumulative total, it needs to be justified by 
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an induction argument. Note that Nn{Q) — Nd(0) = 1, so the values arc correct for 
the first cycle. Since the value is 1 at the end of the first cycle, the pattern repeats 
at each subsequent cycle. Note that we know a priori that < A^jv(i) — ^^{t) < 3 
because the boundary of SG consists of three points,, and this is consistent with 
the values in the Table. 

Note that the multiplicities are the same for the first seven eigenvalues of each 
cycle, but the last three entries in the Table have multiplicities that depend on the 
value of j in (2.1). 

We can also say explicitly what the values of the values of the primitive eigen- 
values are. Define 



„^ / X 5 ± \/25 - 4t 

(2.2) ip±{t) = ^2 and 

(2.3) V(i) = | lim 5VL'=)(f) 

where (/?^'' denotes if- composed with itself k times (the convergence of the limit 
in (2.3) follows from the Taylor expansion (p-{t) = + 0{t^) near t = Q). More 
generally, \{ 5 = {5i,52, ■ ■ ■ ,5n) is a finite sequence of ± signs, we write ips = 
'fSi o </'(52 o • • • o ips^ and \5\ = m. Then 
(2.4) 

{A(f)} = {5V'(p), 52v.(<p+(p)), 5V(<P++(P)), ^^^{<f+-m. 5l^l+V(¥'5(p)), . . .} 

for p = 2,3 or 5, and Si = + but otherwise 6 is unrestricted. Therefore there are 
2m- 1 ciioices of 6 with \6\ = m, and since 5- is increasing and 5+ is decreasing 
it is necessary to use an unconventional ordering in order for the list (2.4) to be 
in increasing order. However we always have 5\^\+^ip{ips{p)) < 5\^'\+^ip{ips>{p')) if 
1^1 < \6'\. This follows from the estimate 

(2.5) ^(3)5l^l+^ < 5l^l+V(¥'5(p)) < ^i'(5)5l^l+^ 

which follows from the estimate 3 < (p+{t) < 5 on < t < 6 and the fact that ip is 
increasing. 

The following is the key Lemma in the proof of the Theorem. Note that it is 
really (1.9) for the special values t = Xn^. 
Lemma: For any n and j, 

(2.6) 7V(5^+2A(f)) =3^'7V(52A(f)). 

Proof. For n = 1 it is known that A''(5-'+^Ap^) = 3-^+^, the total number of vertices 

in the level j + 2 graph approximation Tj+2 of SG, It is easy to see in any case from 
the Table that N{b^x'-^^) = 27 and N{5^X^^^) = 81. We will show by induction on 
k that 

(2.7) N{5^+^xf) = 3N{5^+^xf) 

and this verifies the k = 1 case. So we assume that the result is true for fc — 1. 

The argument is slightly different depending on the parity of fc, so assume first that 
fc = 2£ — 1 is odd. Then fc — 1 = 2-' (2£' — 1) is even, and the induction hypothesis 
is 

(2.8) N{5^'+^xf) = 3N{5^'+^xf). 



EXACT SPECTRAL ASYMPTOTICS ON THE SIERPINSKI GASKET 



5 



What we need to show is 

(2.9) iV(53Af ^) - NiB^'+'^xP) = 3 (^N{5''xf ^) - iV(5-''+2A^?')) . 

Note that on the right side of (2.9) we are counting the last eigenvalue in cycle 
Cfc-i (multiplicity 3-' "'"^ + 1) and the first nine eigenvalues of cycle Ck (multiplicities 

1 + 2 + 2 + 3 + 2 + 2 + 1 + 4 + 9 = 26), so the right side of (2.9) is 3^'+^ + 3*. On 
the left side of (2.9) wo arc counting the last eigenvalue in cycle C2fe-2 (multiplicity 
3^ +2 _|_ oil ^jjg eigenvalues in cycle C2k-i (multiplicities 1 + 2 + 2 + 3 + 2 + 

2 + 1 + 4 + 9 + 4 = 30) and the first nine eigenvalues in cycle C2k (multiplicities 
1 + 2 + 2 + 3 + 2 + 2 + 1 + 10 + 27= 50) for a total of 81 + 3^'+2. This verifies (2.9). 

Finally, assume k = 2^ {2£ — 1) is even, and fc — 1 = 2^' — 1 is odd. The induction 
hypothesis is 

(2.10) N{5^X'^P) = 3N{5^xP), and we need to show 

(2.11) N{5^+^xf) - iV(53A^?^) = 3 (^N{5^+^xf - N{5''xf)^ . 

Wc arc counting the same eigenvalues as before but the multiplicities are different. 
On the right side the multiplicities are 4 and 1 + 2 + 2 + 3 + 2 + 2 + 1 + {3^+^ + 
1) + 3J+2 = 14 + 3-'+^ + 3J+2. On the left side the multipUcities are 10 from 
cycle C2k-2, 1 + 2 + 2 + 3 + 2 + 2 + 1 + 4 + 9 + 4 = 30 from cycle C2k-i, and 
1 + 2 + 2 + 3 + 2 + 2 + 1 + (3^+2 + 1) + 3^+3 = 14 + 3^+^ + 3^+^ from cycle Cafe, for 
a grand total of 54 + 3^+^ + 3^+3. □ 

There appears to be a lot of arithmetic in this proof with no clear explanation. 
Certainly it would be desirable to have a more conceptual proof. 

3. Proof of Theorem 
We consider two intervals between 5'^X^f^ and its neighboring eigenvalues, namely 



(3.1) 




c\(5) c2x(3) 
c:2x(3) ^>(5)^ 



) 

The function N{t) is constant on these intervals, 



(3.2) N{t) 



N{5^xf), t e A'g. 



More generally, N{t) is constant on the intervals 5^A£ and 5^A'g for j > 0, namely 

This uses the Lemma and the fact that the multiplicity of 5^+^A^^^ is This 
gives us (1.9) with mo{t) = on U with 



(3.4) Git) 



i (iV(52AfV9)t-", teAi 
i7V(52Af )t-«, teA',. 
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For any n G Z, wc similarly obtain (1.9) with mo{t) = max{— n,0} for t G 5^{Ae U 
Air ) with G extended by multiplicative periodicity. We also observe from the Table 
that (1.10) holds on the same intervals with 

(3.5) Gi{t) -- 

Now let A = U^-oo Ufci 5" (^£ U A^). To complete the proof we have to 
show that the complement of A has measure zero. Wc can restrict attention to the 

interval [0, 5] because 5A = A. Write A'l = AeU{b^xf^}UA'f = {^^^^^^n-i ' ^A^^^) 

and A" = U^-c» U^li 5" A". The complement of A" differs from the complement 
of j4 by a countable set of points, and it simplifies matters to deal with A". Now 
each of the intervals A'^ has the form 

(5'"V(¥'5(5)),5'"V(^^'(5)) 

where = |(5'| = m — 1 and we no longer require that Si and 6'i be + (this uses 
the fact that tpiip-ix)) = 57p{ip{x)). Then 5-™A^' = (V'(<^5(5), V(w(5)) Q [0,5]. 
The composition with ^ is just a diffcrcntiablc distortion factor, so we drop it and 
define B£ = {(fs{5),(fis'{^)) and B = U^i-^^- problem is to show that the 

complement of B in [0, 5] has measure zero. 

Call this complement J. We claim that J is just the Cantor set defined by 
the mappings (p+ and (p- on [0,5], or equivalently the Julia set for the mapping 
z — )• z{5 — z) since ip± are the the inverses of this mapping. It is known that this 
Julia set has measure zero, but in any case this is easily seen from estimates on the 
derivatives of ip± . 

Note that = (^^,^^) so the complement of Bi is <^_([0, 5])U</j+([0, 5]). 

Then B2 = ifi-{Bi) and B3 = (f^{Bi) so the complement of Bi U B2 U B3 is 

(p ([0,5]) U (^-+([0,5]) U ^++([0,5]) U v3+_([0,5]). Iterating this reasoning, the 

complement of U^=r^ -^^ U|5|=m V'sii'^^ 5]). This verifies the Cantor set descrip- 
tion of J. □ 
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